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Low-Thrust Systems Design 



• Low-thrust trajectory & s/c hardware system are tightly coupled 

- Definition of traj. dependent on propulsion system, LV 

- SEP has variable power & dependent on array size 

• Systems design problem 

- Different possible l sp , power levels, number of thrusters, launch vehicle 

- Realistic engine, array models are discrete 

- Hybrid optimal control problem 

- Design space is multimodal, mixed parameter, often expansive 

• Traditional approaches to sample trade space 

- Directly vary power & l sp in optimization formulation 

- Simplified models, characteristic solutions 

- Parametric studies, grid searches 

• Limitations 

- Trajectory opt. requires initial guess; locally optimal only 

- Only single-objective opt. strategies employed 

- Grid searches intractable 

- Limited fidelity w/out trading realistic hardware models 

- No full mapping of optimal trade space 


BOL p 0 : 30 kW, 3 thrusters 



BOL p 0 : 68 kW, 7 thrusters 

departure 
Bennu 
1/2/2026 
m = 14938 kg 



intercept 
Earth 
4/16/2028 
w* = 1.414 km/s 
DEC = -10.4 
m = 12614 kg 



launch 


Earth 


8/21/2023 


C3 = 1.575 km 2 j s 2 


DLA = 28.5 


m = 10526 kg 



LT rendezvous 
Bennu 
10/29/2024 
m = 8959 kg 




Objective 


IV>SA 


Solve multi-objective, low-thrust systems optimization 
problem to fully map optimal systems trade spade 


Method should be: 

• Capable of global trajectory & systems 
parameter search 

• Automated 

• Free of user-defined initial guess 

• Able to search broad design space 

• Medium fidelity for preliminary design purposes 

• Efficient 
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Multi-objective Optimization 
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• Want to optimize any number of mission design metrics 

- e.g., payload mass, TOF, array size, ref. power, number of thrusters 

- Often coupled & competing 

- Fully map mission trade-offs between optimal solutions 


• Optimize multiple objectives 
simultaneously 

- Entire set of optimal solutions 

- Goal: generate representation of 
Pareto front 

- Traditionally use repetitions of 
single objective technique 
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Multi-objective Systems Optimization 
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Approach: Solve coupled problem simultaneously w/ hybrid optimal 
control algorithm 

• Multi-objective genetic algorithm (GA) as outer loop systems optimizer 
around direct-method inner loop trajectory optimizer 

- Non-dominated Sorting Genetic Algorithm II (NSGA-II) searches over systems 
parameters, defining trajectory problem 

- Monotonic basin hopping (MBH) + sequential quadratic programming (SQP) solves 
trajectory problem 


Initial generation 


Payload 

Mass 



Final generation 


Population evolves via 
genetic operators 


Payload 

Mass 


• • • 


BOL power 


BOL power 
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Genetic Algorithm 


Initial random population 
generated 




Models Darwinian evolution 
- Mimic natural selection & reproduction 
Searches with population of designs 
Globally search design space 
No initial guess required 
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NSGA-II 




• Develops globally-optimal Pareto solutions using non-dominated sorting 

- Evolves population towards Pareto front 

• Fitness assignment based on “nearness” to Pareto front 


- x 1 dominates x 2 if: 

P = 12,-,n 0 bj 

and 

3 P : f p (*! ) < f p (x 2 ) P = 1 ,2,..., n obj 

• If neither design dominates other, they 
are non-dominant 

• Non-dominated sorting : 

- Assign fitness based on design’s non- 
dominated front 

- Designs closer to Pareto front -> better 
fitness & more mating opportunities 
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Low-Thrust Trajectory Optimization 


Need automated, robust method that does not require initial guess 


Solution: apply a global-local hybrid algorithm 

- Formulate problem based on Sims & Flanagan transcription 

- Monotonic basin hopping (MBH) drives global search 

- Gradient-based optimizer solves NLP (SNOPT used) 

Robust & efficient formulation 


Continuous thrust approximated 

- Trajectory discretized into segments 

- Impulsive AV at segment midpoint 

Efficient constraint handling 

- Gradients guide search 

- Robust & efficient formulation 

Proven approach in EMTG software 
(Evolutionary Mission Trajectory Generator) 


O Control Node 
A A Match Point 
^ » Impulse 
I Segment Boundary 



Leg 


From Sims and Flanagan 
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Monotonic Basin Hopping + SQP 



• Stochastic, global search scheme 

• No initial guess required 

• Adept at multi-modal problems w/ clustered local minima 

• Stochastic “hops” evaluated from base solution 

- Pareto distribution balances exploration & exploitation 
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Combine P t and Q t to form R t 


Non-dominated sort to generate P, + 1 

I 

Crowded tournament selection on Pn i 
to generate parent pool S 


Uniform crossover of S to create 
offspring population. Q M 


Mutate Q U1 











Synergistic relationship between outer & 
inner loops 

Generates globally optimal Pareto solutions 
for mission trade evaluation 

Any number of objectives viable 

Flexible to any unique mission constraints, 
trajectory constraints enforced in EMTG 


10 











Conclusi 

ons 





' '■ - 


• Hybrid optimal control algorithm developed for low-thrust spacecraft systems 
design 

- Outer loop: NSGA-II solves systems optimization problem 

- Inner loop: MBH+SQP solves trajectory optimzation 

• Generates globally optimal Pareto solutions for mission trade evaluation 

• Automated 

• Any number of objectives viable 

• Ability to trade discrete, realistic hardware models 

• General applicability to any interplanetary, low-thrust mission 

- Flexible to any unique mission constraints, trajectory constraints enforced in EMTG 

• Can make large systems problems computationally tractable 
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Asteroid Robotic Retrieval Mission: return asteroid boulder or entire asteroid 

- Extensibility option is to return boulder from Deimos 

- Want to understand how return mass & TOF are affected by array size, # of thrusters 

-> Multiple objectives: maximize return mass, minimize TOF, minimize BOL power, minimize # 
of thrusters (all coupled) 


Mission Parameters 


System Design Variables 


Design Variable 

Integer 

Value 

Resolution 

Launch option 

[0,1] 

{Delta IV-H from LV curve, 
Delta IV-H with LG A] 

- 

Solar array size 

[0, 20] 

[30, 70] kW 

2 kW 

Launch window 
open epoch 

[0,4] 

{2020, ...,2029} 

1 year 

Flight time 

[0, 26] 

[700, 3300] days 

1 00 days 

Engine type 

[0,2] 

{high-lsp, medium-thrust, high- 
thrust} 

- 

Number of 
engines 

[0,5] 

[2, 7] 

1 


Description 

Value 

Launch window 

1 year 

Wait time at Bennu 

[430, 700] days 

Min. spacecraft mass with 2 
thrusters 

5991 kg 

Additional dry mass per extra 
thruster 

75 kg 

Max. depart, mass if lunar gravity 
assist (C 3 < 2.0 km 2 /s 2 ) 

11191 kg 

Max. departure mass if direct 
launch (C 3 = 0.0 km 2 /s 2 ) 

10796 kg 

Maximum C 3 if direct launch 

6 km 2 /s 2 

Post-mission AV, l SD 

75 m/s, 3000 s 

Thruster duty cycle 

90% 

Solar array modeling 

1/r 2 

Spacecraft bus power 

2 kW 

Propellant margin 

6% 
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Pareto-Optimal Solutions 
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BOL power at 1 AU (kW) 




Final journey mass increment (for maximizing sample return 


Optimal Trade Space 






6.4 
6.0 

5.6 

5.2 
4.8 

4.4 
4.0 

3.6 

3.2 


• Sharp increase in maximum return mass w/ increasing power 

- Increase in dry mass for increased power not accounted 
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Flight time (years) 






Optimal Design Parameters 
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• Distinct grouping of engine modes based on TOF 

- Return mass plateaus for different engines 
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Backup 
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• Asteroid Robotic Retrieval Mission (ARRM) Option B target 


Mission Objective 

Return a large boulder from Bennu 

Launch Vehicle 



Delta IV Heavy direct (C3 < 6.0) 


Delta IV Heavy with lunar flyby (C3 2.0) 

Power System 


Array power at 1 AU 

chosen by optimizer 

Cell performance model 

1/r 2 

Spacecraft bus power 

2.0 kW 

Power margin 

0% 

Propulsion System 


Thruster 

chosen by optimizer (high-lsp , medium thrust, or high-thrust versions of a large Hall thruster) 

Number of thrusters 

chosen by optimizer (2, 3, 4, 5, 6, 7); dry mass increases by 75 kg for each addtl thruster 

Duty cycle 

90% 

Propellant tank 

unconstrained 

Mission Sequence 

Direct travel to Bennu followed by direct return to C3 2.0 for lunar flyby capture 

Inner-Loop Objective Function 

Maximize sample return mass 

Outer-Loop Objective Functions 

Sample return mass 


Solar array size 


Number of thrusters 


Flight time 
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Bennu Sample Return: Outer-Loop Menu 



Launch Year 

Code 

Year 

0 

2019 

1 

2020 

2 

2021 

3 

2022 

4 

2023 


Power Supply at 1 AU 

Code 

Array Output 

0 

30 

1 

32 

2 

34 

3 

36 

20 

70 


Earth Departure Type 
Code Type 

0 Delta IV-H direct 

1 Delta IV-H w/ LGA 


Flight Time Upper Bound 


Thruster Type 

Code 

Days 

Code 

Thruster 

0 

700 

0 

13 kW Hall (High-lsp) 

1 

800 

1 

13 kW Hall (medium-thrust) 

2 

900 

2 

13 kW Hall (High-thrust) 

3 

1000 



4 

1100 



5 

1200 



7 

1300 



8 

1400 



9 

1500 



10 

1600 



26 

3300 




Number of Thrusters 

Code 

# Thrusters 

0 

2 

1 

3 

2 

4 

3 

5 

4 

6 

5 

7 


102,060 possible 
combinations 
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Final journey mass increment (for maximizing sample return) 


Bennu Sample Return: 

Final Generation Trade Space 
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4.0 
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2.0 
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Number of thrusters 




Final journey mass increment (for maximizing sample return) Final journey mass increment (for maximizing sample return) 


Bennu Sample Return: 


Evolution of Population 
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BOL power at 1 AU (kW) BOL power at 1 AU (kw) 



Final journey mass increment (for maximizing sample ret 


Bennu Sample Return: 
Objective Space 
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Return mass vs. array size 
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Number of thrusters 




Bennu Sample Return: 
Objective Space 


TOF vs. array size 
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TOF vs. number of thrusters 
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Final iournev mass increment (for maximizina sample re 





Boulder Return Mass (kg) 


Bennu Sample Return: 
Optimal Desian Variables 


Departure Type 
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Engine Type 
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Bennu Sample Return: 
Two Trajectories 


A 8-year mission with a 58 kW solar 
array returns a 20 ton boulder 






A 3.3-year mission with a 70 kW solar 
array returns a 2.2 ton boulder 

ie8 „ 
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